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We show that, at sufficiently large strength of the long-range Coulomb interaction, a mass term 
breaking parity (so-called Haldane mass) is dynamically generated in the many-body theory of Dirac 
fermions describing the graphene layer. While the tendency towards chiral symmetry breaking is 
stronger than for the dynamical breakdown of parity at spatial dimension D > 2, we find that the 
situation is reversed at D = 2. The need to regularize the many-body theory in a gauge-invariant 
^ , manner (taking the limit D = 2 — e) is what leads to the dominance of the parity-breaking pattern 

in graphene. We compute the critical coupling for the generation of a parity-breaking mass from the 
' finite radius of convergence of the ladder series supplemented with electron self-energy corrections, 

finding a value quite close to the effective interaction strength for graphene in vacuum after including 
\& . Fermi velocity renormalization and static RPA screening of the Coulomb interaction. 
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I. INTRODUCTION 

During the last years, a great effort has been devoted to the investigation of the so-called graphene, the material 
made of a one-atom-thick carbon layer [lj. Many of the interesting features of graphene arise from the unconventional 
band structure of the carbon sheet, where electron quasiparticles behave like massless Dirac fermions @, |J. The 
relativistic-like invariance of the low-energy electron system is at the origin of outstanding transport properties [J , 
like the reduced influence of scatterers with size above the C-C distance[5| or the perfect transmission through one- 
dimensional potential barriers @. 

In the graphene layer, the Coulomb repulsion between electrons constitutes the dominant interaction. This makes 
the electron system to be at low energies a variant of Quantum Electrodynamics, but placed here in the strong 
coupling regime as the ratio of e 2 to the Fermi velocity vf of the electrons is nominally larger than one. There have 
been already several proposals to observe unconventional signatures of the interacting electrons in graphene, including 
the anomalous screening of impurities carrying a sufficiently large charge It has been also found that the own 

e-e interaction in the layer should lead to a linear dependence on energy of the quasiparticle decay rate fill]. as a 
consequence of the vanishing density of states at the charge neutrality point, and in agreement with measurements 
carried out in graphite [l2j|. 

Yet the effects of electron correlations have been quite elusive in graphene (setting aside the observation of the 
fractional quantum Hall effect). In this respect, a number of theoretical works have studied the dynamical breakdown 
of the chiral symmetry of the Dirac fermions |13l - [27| . being remarkable that most part of the approaches have asserted 
its viability, with a critical coupling below the nominal interaction strength of graphene in vacuum (e 2 /Attvf ~ 2.2). 
Experimentally, it has been observed instead the increase of the Fermi velocity when approaching the charge neutrality 
y—( ' point I28| . reflecting the renormalization of the effective interaction strength towards a vanishing value in the low-energy 
(Nl! theoryjSa H^l . 

It has been proposed that this scaling of the electron system approaching the noninteracting limit could explain the 
absence of any effect of dynamical gap generation in graphene. Some theoretical studies have indeed incorporated the 
electron self-energy corrections in the many-body theory to show that they can push the critical coupling for exciton 
condensation above the largest possible value attained for graphene in vacuum]!?], l31| . Anyhow, the theoretical 
analyses have dealt mainly with the static RPA screening of the Coulomb potential, that tends to underestimate the 
effective strength of the interaction at low energies. The consideration of a more sensible dynamical screening has 
shown to lead to smaller values of the critical coupling for the excitonic instability, below the value e 2 /Airvp = 2.2 

These preceding investigations have focused on the dynamical generation of a gap from the formation of a staggered 
charge density in the graphene lattice, while the possible condensation of other order parameters has been mostly 
overlooked (see however Refs. [H, [33j]). Nevertheless, already in the early studies of Quantum Electrodynamics at 
spatial dimension D = 2, it was realized that the dynamical symmetry breaking could develop in two different ways. 
That is, one could have the condensation of a mass breaking the chiral symmetry but preserving the invariance under 
parity of the Dirac fermions [34l - [3(| . or otherwise a mass maintaining the chiral symmetry but breaking the invariance 
under parity [U l38| . The prevalence of one pattern of symmetry breaking over the other was shown to depend in 
general on the balance between different local four-fermion interactions in the electron system [38|. 
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The two different order parameters characterizing the dynamical symmetry breaking, the parity-invariant mass 
and the parity-breaking mass, are not related by any fundamental symmetry in the graphene electron system, which 
means that they can be generated quite independently. The possibility of introducing a mass operator breaking the 
invariance under parity was considered in the context of a tight-binding model of electrons in the honeycomb lattice 
in the seminal work of Ref . [3^ . The dynamical generation of the parity-breaking mass establishes in fact a bridge 
between graphene and the so-called topological insulators, as one of the consequences of the condensation of the 
parity-breaking order parameter would be the spontaneous development of loop currents in the honeycomb lattice. 

In this work we study the dynamical generation of the parity-breaking mass term in the theory of Dirac fermions 
in graphene, taking the long-range Coulomb repulsion as the relevant interaction between electrons. We will consider 
the dynamical symmetry breaking analyzing the corresponding order parameter in the many-body theory, comparing 
it with the more conventional effect of exciton condensation breaking chiral symmetry but preserving parity. We 
recall that the space of Dirac fermions at D = 2 is defined by the algebra of matrices {~f a } (c = 0, 1, 2) satisfying the 
anticommutation relations 

{ 7 ^ 7 4=2diag(l,-l,-l) (1) 

In this notation, the order parameter for the parity-invariant mass will correspond to the fermion density averaged 
with the 70 matrix, while the generation of the parity-breaking mass will be given by a nonvanishing average of the 
fermion current iesijli'Yj- 

One can check that, at D — 2, the product of all the 7 matrices, 75 = 7 o 7 i 7 2: commutes with any other matrix 
in the space of the Dirac fermions, so that its action is equivalent to the identity. This means that, computing 
for instance the ladder series of many-body corrections, all contributions to the 70 vertex should be formally the 
same as those to the vertex for the current i^uj^fi'Jj- The sum of the ladder series provides a particularly sensible 
approach to dynamical symmetry breaking, as it encodes the most divergent diagrams at each level of the perturbative 
expansion [40]. It is actually the need to devise a consistent regularization of those divergences which plays a central 
role in the evaluation of the different order parameters. 

We will show that, in the theory regularized with a method chosen specifically to preserve its gauge invariance (such 
as dimensional regularization), the perfect match between the order parameters for the parity-breaking and the parity- 
invariant mass is lost. This situation is similar to the well-known case of conventional Quantum Electrodynamics at 
spatial dimension D = 3, where there is no regularization that can simultaneously preserve the gauge and the chiral 
symmetries, which makes unavoidable to discard the latter in favor of the first. 

In the case of graphene, the use of dimensional regularization has shown to be a most suitable choice to preserve the 
underlying gauge invariance of the many-body theory [2^, [U |42[. Recently, it has been checked that such a method 
is singled out as the right approach to reproduce the corrections to the DC conductivity calculated from the lattice 
model [43]. In our problem, the need to regularize the many-body theory while preserving gauge invariance is also 
what produces the anomalous mismatch between the respective vertices for 70 and the current itzijl-Hj- 

We will see that, at D greater than 2, the divergent corrections to the vertex for the parity-invariant mass are in 
general greater that those to the vertex for the parity-breaking mass, while the situation is reversed for D < 2. The 
diagrams in the Dirac many-body theory are made convergent by analytic continuation to spatial dimension D = 2 — e, 
which explains that the tendency to develop a nonvanishing order parameter for the parity-breaking mass becomes 
dominant in the physical renormalized theory at D — > 2. As a consequence, we will see that there is a critical coupling 
for the dynamical breakdown of parity in graphene which turns out to be smaller than the critical coupling obtained 
in the more conventional mechanism relying on chiral symmetry breaking. 

In the next section, we set up the framework to describe the many-body theory of Dirac fermions, highlighting 
the classical scale invariance that opens the way to apply the renormalization group approach. We focus on the 
computation of the vertices for 70 and the current itzij^Cij in Sec. Ill, where we show that it is possible to define 
cutoff-independent observable quantities from the sum of ladder diagrams. We find that the two ladder series for the 
renormalized vertices have finite radii of convergence, that we compute in each case to obtain the respective critical 
couplings signaling the dynamical generation of the parity-invariant and the parity-breaking mass. Finally, we discuss 
in Sec. IV the physical implications of the dynamical symmetry breaking arising from a parity-breaking mass. 

II. DIRAC MANY-BODY THEORY 

The low-energy electron quasiparticles are disposed in graphene into conical conduction and valence bands that 
touch at the six corners of the Brillouin zone|4J. There are two inequivalent classes of electronic states, that can be 
disposed around a pair of independent Fermi points. Thus, the low-energy electronic excitations can be encoded into a 
couple of four-dimensional Dirac spinors {ipi}, which are characterized by having linear energy-momentum dispersion 
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e(p) = Uf|p|- The index i accounts here for the two independent spin projections, but it will be omitted to simplify 
the notation in what follows. The kinetic term of the hamiltonian in this low-energy theory is given by 

H = -iv F J d 2 r ^(r) 7 • W(r) (2) 

where ip = ^'70 and {7,7} is a collection of four-dimensional matrices satisfying the relations ((T|). It is convenient to 
represent them in terms of Pauli matrices as 

70,1,2 = (o"3, 03 <7i, 0-30-2) ® r 3 (3) 

where the first factor acts on the two sublattice components of the graphene honeycomb lattice and the second factor 
operates on the set of two independent Fermi points. 

The dominant e-e interaction is given in graphene by the long-range Coulomb repulsion. This interaction does not 
have a conventional screening at long distances, given the vanishing density of states at the Fermi points connecting 
the conduction and valence bands. A suitable description of the Coulomb interaction starts then by considering the 
unscreened potential V(r) = e 2 /47r|r|. The full hamiltonian including the e-e interaction can be represented as 

H = -iv F [ d 2 r^(r) 7 • V^(r) + |- / d 2 r x [ d 2 r 2 p(n)- -p(r 2 ) (4) 

J 8tJ J |ri — r 2 | 

with p(r) = V^rhoVKr). 

The prevalence of the Coulomb interaction at low energies can be understood from the scaling properties of the 
system governed by ((4]). The long-range Coulomb repulsion is the only interaction that is not suppressed, at the 
classical level, when scaling the many-body theory towards the limit of very large distances or very low energies. 
From Q , the total action of the system is given by 



S = j dt j d 2 r Tjj(r)(ij d t + iv F j • V)ip(r) 

This action is invariant under the combined transformation of the space and time variables and the scale of the fields 

t' = st , r' = sr , V' = (6) 

This means in particular that the strength of the Coulomb interaction is not diminished when zooming into the 
low-energy limit s — > 00. One can check that any other e-e interaction without the l/|r| tail, as those that arise 
effectively from phonon exchange, would be suppressed at least by a power of 1/s under the change of variables ([5]), 
implying its irrelevance in the low-energy limit. 

On the other hand, the many-body theory does not preserve in general the scale invariance of the classical action 
((5|) , as a high-energy cutoff has to be introduced to obtain finite results in the computation of many-body corrections 
to different observables. The analysis of the cutoff dependence of the many-body theory provides deeper insight into 
the effective low-energy theory. If the theory is renormalizable, it must be possible to absorb all powers of the cutoff 
dependence into a redefinition of the parameters in the action ([5]) . This should be therefore modified to read 



S = Z kin J dt J d 2 r ^(r)(z 7o 9 t + iZ v v F j • V)V>(r) 

e 2 r r r i 



-Z int — J dt J d 2 n J d l r 2 P ( ri ) p(r 2 ) (7) 

The assumption is that Z^i n ,Z v and Z mt (and other renormalization factors for composite operators that do not 
appear in ([5])) can only depend on the cutoff, while they must be precisely chosen to render all electronic correlators 
cutoff independent [44J • 

This possibility of absorbing all the divergent dependences on the cutoff in a finite number of renormalization factors 
has been checked in the above theory in the limit of a large number N of Dirac fermion flavors [30ll45j . as well as in the 
ladder approximation to vertex diagrams supplemented by electron self-energy corrections [42j. We will use this latter 
approach in the computation of the vertices for 70 and the current ie3ijjijj ■ We recall that in the ladder approximation 
there is no need to renormalize the electron quasiparticle weight, so that Z^ = 1, while the renormalization factor 
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Z v for the Fermi velocity coincides with the expression obtained at the first perturbative level 42]. The computation 
of the expectation value of other fermion currents requires however the introduction of new renormalization factors 
in order to obtain finite, cutoff- independent results, as shown in Ref. [42| . These renormalization factors encode in 
general valuable information about the possible condensation of different excitations and the corresponding dynamical 
symmetry breaking. 

The gapless character of the electronic spectrum is in principle protected by two different symmetries of the hamil- 
tonian Q. This is invariant under a Zi chiral symmetry (CS) that consists in the exchange of the two-dimensional 
spinors attached to the two independent Fermi points 

C75:^(r)-»-(l®7i)V(r) ( 8 ) 
We could add to the hamiltonian (j4|) a term breaking the invariance under ([8]) 

H rn = m J d 2 r ^(r)^(r) (9) 

while respecting the space-time symmetries (46|. This term, corresponding to a conventional parity- invariant mass, 
would certainly open a gap in the electronic spectrum. But this same effect can be obtained with a different pertur- 
bation preserving the chiral symmetry (JSJ 



H m > = ml J d 2 r ip(r) eujlili ip(r) (10) 

One can check that (fit))) breaks instead the invariance under parity (P) 

P :ip(x,y) -> (crt®Ti)ip(x,-y) (11) 

The term (fTU|) corresponds to the parity-breaking mass mentioned before [39]. It can be shown that there are no other 
operators (not implying the electron spin or the exchange of excitations about the two Fermi points) that may open 
a gap in the electronic spectrumJU Eg. We finally recall that, while ((9]) or (flQ|) may not be present nominally in the 
hamiltonian, they can be generated in the many-body theory at sufficiently large interaction strength, leading to the 
dynamical breakdown of the chiral symmetry or parity in each respective case. 



III. RENORMALIZATION OF STAGGERED CHARGE DENSITY AND LOOP CURRENT 

OPERATORS 



We focus then on the possible development of a vacuum expectation value for the parity invariant operator 

p 3 (r)=?(r)V(r) (12) 

as well as for the parity-breaking bilinear 

p 3 (r)=?(r) 7o lV/Hr) (13) 

with 

^3 = iesijUli (14) 
We will characterize such a phenomenon by inspection of the respective vertices 

r 3 (q,w, z ;k,wfe) = (p 3 (q, uj q )ip(k + q, uj k + w 9 )^(k, Wfc))iPi (15) 

and 

r 3 (q,u; 9 ;k,u; fe ) = (p 3 (q, uJ q )4>(k + q, uj k + ui q )ip(k, u>fc))ipi (16) 

where 1PI denotes that we take the one-particle irreducible part of the correlator. We are going to see that the above 
vertices have a finite radius of convergence in perturbation theory, which is the signature of the dynamical breakdown 
of symmetry driven by the interaction. 

A sensible approach to determine the critical interaction strengths at which (|15j) and (|16|) diverge consists in 
performing the sum of ladder diagrams for the two vertices. The ladder contributions contain indeed the most 
divergent part of the vertex at each level of the perturbative expansion 40] . The ladder series built in that way can 
be easily encoded in the self-consistent diagrammatic equation represented in Fig. [1] In what follows, we will see how 
to obtain the respective critical couplings in a manner that they do not depend on the high-energy cutoff needed to 
regularize the theory, relying only on observable quantities of the electron system. 
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FIG. 1: Self-consistent diagrammatic equation for a generic vertex r, in the ladder approximation. 

A. Loop current vertex 

We study first the vertex for the loop current operator. In order to make the approach covariant and valid in any 
spatial dimension, we consider the correlations of the loop operator 

Pij( T ) = ^(r)7o*[7i,7i]V>(r) (17) 
This is now a second-rank tensor, leading to the vertex 

T i j(q,u q ;'k.,u)k) = (pij(q,0J g )ip(k + q,oj k + w 9 )-0(k, Lo k ))iPi (18) 

In the ladder approximation, the vertex rjj(q, w g ;k, u>k) has to satisfy in the limit q = and u q — the self- 
consistent equation 

*7o 7t>7j +»T / TTTd^ - 7o 2 , 2 2 -Iy 0, 0; p, cj p ) 2 -7 0]1 , (19) 

2 J (27r)^ 27r — + w^p^ — z?7 — w~ + i>£p^ — 17/ |k — p| 

where we have regularized the momentum integrals by continuing the number of spatial dimensions to D = 2 — e. 
Accordingly, eo is now a dimensionful coupling which demands the introduction of an auxiliary momentum scale fi 
through the relation 

e = fi e/2 e (20) 

At this point we may resort to a perturbative resolution for the vertex, noticing that it can only depend on the 
dimensionless ratio /i/|k|. It is easily realized that the operator with the tensor [74, 7y] can mix with the other 
antisymmetric second-rank tensor structure one can build, corresponding to ^ikj — ki'jj. Thus, the most general 
solution of Eq. (fH)]) can be expressed as a power series in the coupling Ao = Cq/Atyvf 

r iJ (0,0;k, Wfc ) = z 7 ob,7,] (l + E A 0|k^J +Ho( W -^)(7-n)^A^ (21) 

where we have called n = k/|k|. 

Each term in the series (|21l) can be obtained by means of an iterative procedure from the previous orders. By 
inserting a term proportional to 7o[7ii7j] a t the right-hand-side of Eq. (|19p . we get for instance 



r 2 



d D p duJp -ijoujp + v F j p 1 -ijouip + v F j ■ p 1 
7o =2 , .. 2 o 7o [7i - 7j] ^ , , 2 7o 7j— 



(2tt) c 2tt ,u t^ + t^p 2 ,UL ' 4 " JJ |p|- UJ 2 p + v 2 F p 2 

= A 7o[7» i 7j] |k ^ + e i)e + A 7o(7 i n j - n zlj )h ■ n )p^ 



(22) 



with 



_ ^ r(^) r (^)r( ¥ ) 

4 ^(M)r 2 -^ 



2 

(4lr) . /a r(i + ^)r(^)r(^) 
" w 2 V H-(» ¥ .)r(2-l^*' 
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Moreover, when a term proportional to Joili^j — ^i7j) is introduced at the right- hand- side of the integral equation, 
we obtain 

d D p duJ p -ijouJ p + v F j ■ p 1 -i7 aJp + v F j ■ p 1 
"7o ^ 2 I ,. 2 „9 7o(7;™j " ^7i)(7 • nJo^r ^2 , .. 2 To, 



2 7 (2tt) d 2tt ,u U2 +W 2p2 '| p |- S72 + U 2 p2 ' u |k-p| 

1. r , 4(£) 1 B„(e) 
- 2 A °7o[7i,7j] |k|(n+1)6 + 2 A o7o(7i^ - n. i7j )( 7 • n) |k|( „ +1)e 



(25) 



The above formulas can be summarized in the recursion relations 

t n +i 



A„(e) [t n + \t[)j (26) 
4+i = B n {e) (t n + \A (27) 



In compact form, one can also write the relation 

tn+i + \t' n+1 = U n (e) + \B n {e) \ (t n + ht'^j (28) 

We observe that the result of computing the integral in (|19[) always diverges in the limit e — > 0, leading to a sequence 
of higher-order poles in the e parameter as we look at higher perturbative levels in the solution of the self-consistent 
equation. Good news are however that all these divergences can be reabsorbed after a multiplicative renormalization 
by a single factor Zh , ending up with a finite vertex in the limit e — > 

^ij,rcn — ^h^ij (^^) 

The renormalization factor Zh can have in general the structure 

" k(A) 



Zh = l+ J2'^l (30) 



The important point is that it can be shown that Tij^ Ten can be made finite at e = with a set of functions hi(X) that 
do not depend on the external momenta of the vertex. This is the hallmark of renormalizability, by which one can 
reabsorb the cutoff dependences into the renormalization factors for a few local operators of the theory. 

Most interestingly, this procedure of renormalization works in the present case even when self-energy corrections 
are considered in the internal electron and hole states of the vertex. These supplementary contributions to the ladder 
series can be taken systematically into account by replacing the constant vp in Eq. (flUl) by the effective Fermi velocity 
vf(p) found after dressing the Dirac propagator with the electron self-energy correction S 

7ou; — vp~f ■ p — > jqlu — vpj ■ p £ (31) 

Recalling the expression of the electron self-energy in the same ladder approximation (42j, we get 

MP) -vf + ^(4*) ^ (32) 

The effective vp(jp) is in principle singular in the limit e — > 0, but the divergence can be absorbed by a simple 
renormalization of the Fermi velocity 

VF = Z v VF,ren (33) 

Taking the simple pole structure 

Z v = l + bJ- (34) 
one needs to set b\ = — e 2 /I6irvp^ ren in order to render Vp(p) finite at e -> as a function of i>_F ire n- 
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The remarkable point is that the vertex Tij^ cn , supplemented with the electron self-energy corrections, can be made 
also finite in the limit e — > with a modified set of coefficients in Q30p , expressed now as functions of the renormalized 
coupling 

A = (35) 

The first orders in the analytic computation of the residues hi(X), including the effect of the electron self-energy 
corrections, turn out to be 

1 l + 41og(2) -3 + 7r 2 + 121og(2)(3 + 101og(2)) ^ 

klW = "2 A 32 A 1152 A 

9 + 3tt 2 - 361og(2) + 207r 2 log(2) + 504 log 2 (2) + 1376 log 3 (2) + 30C(3) 4 

12288 

' (13tt 4 + 6tt 2 (-3 + 41og(2)(21 -I- 86 log(2))) 



737280 

+6(51 + 162C(3) + 41og(2)(45 + 21og(2)(-45 + 41og(2)(235 + 5371og(2))) + 222C(3)))) A 5 + . . . (36) 
h f\\ 1 , 2 , l + 41og(2) , 3 -6 + 5^ 2 + 121og(2)(21 + 621og(2)) N4 

h2{X) = 16 A + 96 A + 18432 A 

tt 2 (37 + 2201og(2)) + 3(17 + 41og(2)(-7 + 2 log(2)(299 + 7161og(2))) + 90C(3)) 6 



368640 



A 5 + . . . (37) 



l + 41og(2) . 12 + 7 r 2 + 1561og(2)+3601og 2 (2) ^ 5 

^ (A) - 3072 - A 184320 A + "" (38) 

h 4 (X) = _l + 41og(2) A5+ (39) 

y ' 30720 ^ ' 

The function h\ (A) has particular significance since it can be used to determine the anomalous dimension of the 
vertex, which comes from the dependence on the dimensionful parameter /j, of the renormalized theory [441] 

f ij, Te n ~ M 7 " (40) 

The exponent jh can be defined in terms of the renormalization factor Zh by 

lh = z~ h l^ (41) 



When applying Eq. (|4ip. it is a highly nontrivial result that all the pole contributions coming from (f30|) may cancel 
out to give a finite answer for 7^ at e —± 0. This happens whenever one can enforce the set of conditions [49| 

4rh i+ i(X) - hi(X) -^i(A) + 61(A) -jrhi(X) = (42) 
aX d\ dA 

In that case, the anomalous exponent jh is given by 

7h = -A^i(A) (43) 

One can see that the analytic expressions p6p - ([3^|) satisfy the conditions (l4"2"l) . and we have checked that these 
hold also for the numerical solution of the functions hi(X) that we have obtained up to order A 22 in the perturbative 
expansion. The characterization of dynamical symmetry breaking from 7^ turns out to be then very convenient since, 
according to (|43p . the exponent is a scale invariant quantity depending only on the renormalized coupling A. It is 
actually observed in the numerical resolution that the perturbative orders of /11(A) have a geometric growth, as shown 
in Fig. [21 We can then use the same scaling technique as in Ref. [42}], which allowed to reproduce there the value 
of the critical coupling obtained for chiral symmetry breaking from the resolution of the gap equation 23] . From the 
expansion up to order A 22 , we get now the result for the finite radius of convergence of jh 

X c « 0.508 (44) 

One has to keep in mind that this singularity at A c implies the divergence of the correlators of the loop current 
operator P3(r), as they all need to be renormalized by factors of Zh- This implies that their anomalous dimensions arc 
given in general by multiples of jh ■ We find therefore at A c the onset of a new phase of the electron system that must 
correspond to dynamical symmetry breaking, as it is characterized by the appearance of a nonvanishing expectation 
value of the operator P3 (r) at that critical coupling. 
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FIG. 2: Plot of the absolute value of the coefficients hi (upper points) and di (lower points) in the respective expansions of 
fti(A) and di(A) as power series of the renormalized coupling A. 



B. Staggered charge density vertex 

We review the computation of the vertex ([TBI , which has been addressed in Ref . [27j , in order to establish here the 
comparison with the loop current vertex. In the ladder approximation, r 3 is proportional to the unit matrix, and the 
resolution of the equation in Fig. [T] proceeds by noticing that 



-j Up + vfJ ■ P r /n -7o^p + vfI ■ P r 3 (0,0;p,q; p ) 

^° — — t~ 2 — -r 3 (o,0;p,w p — 9 1 2 2 — -70 = 2^ — 22- ( 45 

—UJ~ + U^p^ — 177 — W 2 + Uj,p J — «77 — W 2 + V F p^ — IT] 

We arrive then at the self-consistent equation 



r,(o, 0; k, Uk) = 1 - 4 / (§r$ r 3(°' <>; p. ^)=J^h5P 



(46) 



(47) 



i 2 — ir7 |k — p| 

The solution of Eq. (|46[) does not depend on the frequency ujk , and the self-consistent equation becomes 

r 3 (0,0;k, Wfe ) = l + ^ f £P T (o . pu ) 1 1 

4 J {2tt) 1j v F \p\ |k - p| 

We resort as before to a perturbative resolution for the vertex, expressing the solution of (|4"T|) as a power series in the 
coupling Ao = Cq/Attvp 

00 

r 3 (O,0;k, Wfe ) = l + ^A^ (48) 

n=l 1 1 

Successive terms in f|48[) can be obtained by means of an iterative procedure, using the formula 



where 



We get therefore the relation 



4 J (2tt) d |p|'" v F \p\ |k-p| u |k|(«+i)^ 



(4,)^r(^)r(^)r m 



Sn+1 — Cn( e ) s n (51) 



It is remarkable that the recursion (|5"Tj) almost coincides with ([25]) . the deviation being a term that vanishes in the 
limit e — > 0. We have in fact the result 

Cn(e) = y4(e) + \B n {e) (52) 
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The agreement between the two recursion for r 3 and T^- would be perfect if we had a factor D/2 multiplying the 
right- hand-side of (f2l)|). which is instead missing due to the second-rank character of the tensor renormalized there. 
It can be shown that the precise consideration of this tensor character is a must for the correct implementation of 
a gauge- invariant renormalization, as observed for instance in the case of the electron current vertex 42] . We are 
going to see that the mismatch in the two ladder series, when combined with the effect of the electron self-energy 
corrections, leads to significant differences in the renormalization factors of the vertices T3 and Tij. 

It can be checked that a multiplicative renormalization by a factor Z rn is enough to absorb the divergences in (|48p . 
Thus, we define the renormalized vertex, finite at e — > 0, by 

r3,rcn = Z m T^ (53) 

The pole structure of the renormalization factor takes the form 

Z m = (54) 

i—l 

At this point, we can consider again the effect of the electron self-energy corrections on the vertex by simply 
replacing the constant Vp inside the integral in Eq. (|47p by the effective Fermi velocity (f3"2"j). After renormalizing 
vf by the factor (f3~4"|). it becomes possible to absorb all the poles at e = in the vertex by an appropriate choice of 
coefficients di(X) depending only on the renormalized coupling A. The first terms in the perturbative expansion can 
be computed analytically, with the result that 

= _1 _ 1 , 2 + 1201og 2 (2) _ 10^1og(2) + 6881og3(2) + 15 C (3) 

y ' 2 8 BW 1152 6144 



13tt 4 + 20647r 2 log 2 (2) + 144 (7161og 4 (2) + 371og(2)C(3)) . 



737280 

1^2,1 ,3 , 57r 2 + 7441og 2 (2) 



X + ... (55) 



d 2 {\) = — A + — log(2) A 



16 24 ov ' 18432 
IIOtt 2 log(2) + 8592 log 3 (2) + 135_C(3) 
184320 



A + . . . (56) 



4(A) = -J-l g(2)A^- 7r2 + 36010g2 ^A^ + ... (57) 
iy ' 768 BV ' 184320 v ' 

<kw = -^ lo ^ x5 + --- ^ 

The most important effect of the renormalization is again the anomalous scaling of the vertex r 3 rcn , arising from 
its dependence on the momentum scale /j, [jij . 

IXrcn - M 7m (59) 
The anomalous dimension j m can be obtained from the renormalization factor as 

7m = 7T- (60) 

This equation turns out to have a finite limit e — > provided the set of conditions 

-£-di+i(A) - di(X) -^d!(A) + h(X) Adi(A) = (61) 
aX aX dX 

are fulfilled 49] . In that case, the anomalous exponent is given by 

7m = -A-^-di(A) (62) 
dX 

therefore depending solely on the renormalized coupling A. 

Quite nicely, it can be seen that the analytic expressions in (|55|) - ([58|) satisfy the equations ([BT]) . We have also checked 
by numerical computation that the perturbative expansions of the functions (fj(A) fulfill that set of conditions, at least 
up to the order A 24 we carried out the calculation. We can therefore rely on Eq. (|6"2"|) to obtain a sensible result for 
the anomalous dimension 7 m . The important point is that, as well as in the case of jh, the power series in A turns 
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out to have a finite radius of convergence. This can be already noticed in the representation of the first orders of the 
expansion for rfi(A), plotted in Fig. [5] The accurate fit achieved in Ref. |42| (from the expansion to order A 24 ) allows 
to get the value for the critical coupling 

X' c « 0.5448 (63) 

The critical coupling A' c should mark the onset of a new phase in the electron system, characterized by a nonvanishing 
expectation value of the operator ^(r). We observe however that the critical value (|63]l is larger than the result (|44|) 
obtained for the coupling A c . This implies that the susceptibility to the development of loop currents should be stronger 
than that for the formation of a staggered charge density in graphene. Consequently, the phase with condensation of 
circular currents should take place first, as long as the system has a relative interaction strength above the critical 
coupling given by (14"4"1) . 



IV. DISCUSSION 



We have seen that, at sufficiently large interaction strength, it is possible to dynamically generate a mass term 
breaking parity (the so-called Haldane mass) in the many-body theory of Dirac fermions describing the graphene 
layer. We have shown that the critical coupling for the development of such a mass is below that corresponding to 
the usual excitonic instability, so that the dynamical breakdown of parity should prevail in graphene over the opening 
of a gap in the electronic spectrum from chiral symmetry breaking. 

As already mentioned, the mismatch between the critical couplings for the dynamical breakdown of parity and 
chiral symmetry may be paradoxical, given that the commutation properties of the vertices for 70 and the current 
ie3ij"fzlj are the same for Dirac fermions at D = 2. We have clarified that the difference between the respective many- 
body corrections to the two vertices arises from the need to regularize high-energy contributions in a gauge-invariant 
manner, by analytic continuation of the momentum integrals below dimension D = 2. The suitability of this method 
to produce precise physical results for the carbon layer has been recently checked in Ref. In this regard, the 
different strength of parity and chiral symmetry breaking is a quantum field theory "anomaly" arising as a remnant of 
high-energy contributions to the low-energy effective theory of graphene, similar to the anomalous non-conservation 
of the chiral current in conventional Quantum Electrodynamics at D = 3 when a gauge invariant method is used to 
regularize the theory. 

It is interesting to observe that the above computation of anomalous exponents can be also applied to investigate 
the dynamical symmetry breaking of the many-body theory of Dirac fermions at any spatial dimension D. This 
is so as the classical scale invariance (JSJ can be extended to hold for the action of the Dirac fermions interacting 
with the same l/\r\ Coulomb potential at arbitrary D. It can be shown then that all cutoff dependences can be 
always absorbed into a finite number of renormalization factors, at least in the ladder approximation supplemented 
with electron self-energy corrections. The anomalous exponents of the order parameters corresponding to 70 and the 
current Ksijlilj display in general finite radii of convergence in our ladder approach, which we have computed for 
spatial dimension 1 < D < 3. The results are shown in Fig. [3] We observe for instance that the critical coupling 
for the conventional excitonic instability is the lowest at D = 3, while we find the reverse situation for D < 2. At 
D = 2, it is actually the need to regularize the many-body theory by computing in the limit D = 2 — e which makes 
the dynamical breakdown of parity to prevail over chiral symmetry breaking at the integer physical dimension. 

Regarding the actual value of the critical coupling (14"4"1) . we may ask whether there is any chance to observe the 
effects of a parity-breaking mass in real graphene samples. At this point, we note that the coupling A c cannot be 
compared directly to nominal values of the graphene fine structure constant, as A should correspond in our approach 
to the effective coupling after screening of the interaction. That is, the value of A c must be referred to the relative 
strength of the dressed Coulomb potential. Under the assumption of static RPA screening of the interaction, we get 
for instance the relation between the effective coupling A and the nominal value a of the graphene fine structure 
constant 

= TT^a M 

The above formula takes into account the effect of two fermion flavors with different spin projection. In the case of 
graphene isolated in vacuum, the bare value a ~ 2.2 leads to A ~ 0.49, which is very close to the computed critical 
value (|4"4"]) . We recall anyhow that the use of the static RPA tends to underestimate appreciably the relative strength 
of the dressed interaction. The value of A obtained after dynamical screening of the interaction has been shown to 
be (starting from a « 2.2) above the critical coupling A^, in (pjB"]) [25|, [53] ■ Therefore, we are led to conclude that the 
electronic state of graphene in vacuum should correspond to the phase with dynamical breakdown of parity above A c . 
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FIG. 3: Plot of the critical couplings A c for the dynamical symmetry breaking of parity (points on the dashed line) and for 
the breakdown of chiral symmetry (points on the full line) depending on the spatial dimension D of the interacting electron 
system. 



We note finally that the phase with the dynamical breaking of parity may not have necessarily insulating character, 
as it corresponds to the spontaneous development of an order parameter (a vector) with two possible orientations. 
We may think of situations in which the space is divided into different domains where the order parameter alternates 
between opposite directions. In these configurations, electronic transport may be possible along the interface between 
neighboring domains. The study of the dynamics of this kind of disorder at a microscopic level is beyond the scope 
of this work, but it is plausible that, for suitably small domains, the phase with the dynamical breaking of parity 
may still have metallic properties. The development of an ordered phase over macroscopic regions may require some 
external source contributing to the complete alignment of the order parameter. This can be achieved for instance by 
introducing a transverse magnetic field, setting in this way a preferred orientation in the parity-breaking mechanism. 
Signatures of the onset of the ordered phase may have been already seen in the metal-insulator transition observed in 
graphene at relatively low magnetic fields [13, [H|. Further efforts should be devoted to the detection of the broken- 
symmetry phase, which should anyhow appear in the close proximity (about 1 meV or below) around the charge 
neutrality point. 
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